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Abstract. Recent PHOBOS measurements of the excitation function for 
the pseudo-rapidity dependence of elliptic flow in Au+Au collisions at RHIC, 
have posed a significant theoretical challenge. We show here that these data 
' are described by the Buda-Lund model. A universal scaling curve, predicted 

\Q , by the Buda-Lund hydro model, describes not only PHOBOS data, but also 

recent, detailed PHENIX and STAR elliptic flow data. This is a consequence 
of perfect fluid hydrodynamics as implemented in the Buda-Lund hydro model. 
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Introduction, terminology. One of the unexpected results from experiments at 
the Relativistic Heavy Ion Collider (RHIC) is the relatively strong second harmonic 
' moment of the transverse momentum distribution, referred to as the elliptic flow. 

Measurements of the elliptic flow by the PHENIX, PHOBOS and STAR collabo- 
rations (see refs. El El EJ El) reveal rich details in terms of its dependence 
on particle type, transverse (pt) and longitudinal momentum (rj) variables, and on 
the centrality and the bombarding energy of the collision. In the soft transverse 
momentum region (px < 2GeV/c), the measurements at mid-rapidity are found 
to be reasonably well described by hydrodynamical models [El El- By contrast, 
differential measurement of the pseudo-rapidity dependence of elliptic flow 1)2(1]), 
and its excitation function have resisted several attempts at a description in terms 
of hydrodynamical models (but see their new description by the SPHERIO model [ 
El). Here we show that these data are consistent with the theoretical and ana- 
lytic predictions that are based on Eqs. (1-6) of Ref. [ 9 , that is, on perfect fluid 
hydrodynamics. 

There is a confusion in the literature if perfect or perhaps an ideal fluid has 
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been produced in Au+Au collisions at RHIC. The BNL press release [El uses * ne 
terminology "perfect liquid" . On the other hand, some other authors, for example 
refs. [ 111! IT2*) are using a different terminology: "ideal fluid". 

Although in many cases these different terms refer to the same set of equations, 
let us clarify that in fact there is a distinction between ideal and perfect liquids 
and the corresponding hydrodynamical equations. The expression "ideal fluid" is 
frequently used either f : for an inviscid, incompressible fluid, or 2: for an inviscid 
fluid. These definitions of "ideal fluid" are inequivalent, and one should note that 
definition 1 that includes the incompressibility is more frequent, as it is part of 
the primary physics courses: "an ideal fluid keeps its volume but conforms to the 
outline of its container" . 

Let us also consider the definitions of perfect liquids. A perfect liquid is defined 
also in two different ways: f: A fluid which looks isotropic in its rest frame is 
called a perfect fluid. 2: A fluid which has no shear stresses, viscosity or heat 
conduction is a perfect fluid. It turns out that these two definitions of the perfect 
fluids are equivalent. Perfect fluids are frequently used in general relativity to 
model idealized distributions of matter. This terminology, with clear definition 
of the energy-momentum tensor that is diagonal in the rest frame of the matter, 
has also been adopted by the Particle Data Group [El- When considering these 
questions on terminology, it is important to note that lack of viscosity does not 
yet imply that a fluid is perfect, as the lack of heat conduction and the vanishing 
of shear tensor needs also be discussed. Second, it is clear that we are discussing 
the time evolution of an exploding fireball, so an ideal fluid of a time independent 
volume is irrelevant for this problem, in contrast to def. 1 of ideal fluids. 

In conclusion, we will continue to use the well defined and unambigous termi- 
nology of "perfect fluid" hydrodynamics in context of RHIC physics, and strongly 
recommend it also for other authors to avoid confusion on this subject. 
The Buda-Lund hydro model. The tool used to describe the pseudorapidity- 
dependent elliptic flow is the Buda-Lund hydro model. This hydro model [ 1141 H3] 
is based on exact analytic hydro solutions discussed in Refs. [El El- Buda-Lund 
type of exact, parametric hydro solutions are being explored at the moment as well. 
General, non-relativistic, ellipsoidally expanding solutions [I16|. and relativistic but 
non-accelerating solutions [ 1171 I18| are known. An important missing link is the 
family of finite, relativistic, accelerating solutions which currently being pursued [ 

El 

The Buda-Lund hydro model is successful in describing the BRAHMS, PHENIX, 
PHOBOS and STAR data on identified single particle spectra and the transverse 
mass dependent Bose-Einstein or HBT radii as well as the pseudorapidity distribu- 
tion of charged particles in Au + Au collisions both at </snn = 130 GeV [EOj and at 
^/snn = 200 GeV [12]. The model is defined with the help of its emission function; 
to take into account the effects of long-lived resonances, it utilizes the core-halo 
model [ I22j , and characterizes the system with a hydrodynamically evolving core 
and a halo of the decay products of the long-lived resonances. 

The elliptic flow would be zero in an axially symmetric case, so we developed 
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an ellipsoidal generalization of the model that describes an expanding ellipsoid with 
principal axes X, Y and Z . Their derivatives with respect to proper-time (expansion 
rates) are denoted by X, Y and Z. The generalization goes back to the original 
one, if the transverse directed principal axes of the ellipsoid are equal, ic X = Y 
(and also X = Y). The deviation from axial symmetry can be measured by the 
momentum-space eccentricity, 

= x 2 -y 2 

£p X 2 + Y 2 ' 

The exact analytic solutions of hydrodynamics [OH 1231 124) . which form the basis of 
the Buda-Lund hydro model, develop Hubble-flow for late times, ie X -^> T ^oo Xr, 
so the momentum-space eccentricity e p nearly equals space-time eccentricity e. It 
turns out that the data fitting is easier if instead of X and Y the eccentricity 
variable e is introduced and complemented by the average radial flow defined as 



u t ) = 1/2/ X 2 + 1/Y 2 . Obviously, the expansion rates are given as and X 



(ut)/y/i - e p and Y = (u t )/y/l + e p . 

Let us introduce A rj additionally. It represents the elongation of the source 
expressed in units of space-time rapidity. Let us consider furthermore that at the 
freeze-out r A?y = Z and Z rs Z r, and so Arj ps Z . Hence, in this paper we extract 
space-time eccentricity (e), average transverse flow (u t ) and longitudinal elongation 
(Arj) from the data, instead of X, Y and Z. 

Although the proof of Ref. [|25] was given only for the case of exact parametric 
solutions of non-relativistic hydrodynamics, we expect a similar behavior at the end 
of relativistic expansions just as well, because the acceleration of the scales (X, Y, Z) 
is due to pressure gradients. However, the pressure hence its gradients as well tend 
to vanishing values in the late stages of both the non-relativistic [ E3 1241 l2~5] and 
relativistic [ED Ej exact parametric solutions of hydrodynamics. 

In the time dependent hydrodynamical solutions, these values evolve in time, 
however, it was show in Ref. [ I25| that X, Y and Z, and so e, u t and A77 become 
constants of the motion in the late stages of the expansion. 

Elliptic flow from the model. The result for the elliptic flow is the following 
simple universal scaling law: 

= f^, (2) 
Iq{w) 

where /„ is the modified Bessel function of order n. For other flow coefficients, the 
prediction is also simple and universal: 

V2n = I n (w)/I (w), (3) 
V2n+1 = 0. (4) 

The model thus predicts a scaling law, valid under certain conditions detailed in 
Ref. [I15|. fulfilled by these data : every V2 measurement falls on the same universal 
curve of I1/I0 when plotted against w. 



4 



M. Csanad, T. Csorgo, R. A. Lacey and B. Lorstad 



This means, that Vi depends on any measurable quantity only through the 
scaling parameter w. In practice this implies that the detailed and experimentally 
known dependence of the elliptic flow on transverse or longitudinal momentum, 
center of mass energy, centrality, type of the colliding nucleus etc is only apparent. 

Before giving the definition of the scaling variable w it is educational to in- 
troduce the effective temperatures or slope parameters in the impact parameter 
direction (T x ) and in the direction out of the reaction plane (T y ), as follows: 



T x = T Q + rn t X 2 T l 2 , (5) 

Tn 



T y = T Q +m t Y 2 " 2 , (6) 
and 

rat = m t coshes - y). (7) 



Here m t = ym 2 + p x +p y is the transverse mass, a = (AT/T) t = (T — T s )/T s 

measures the temperature gradient in the transverse direction, at the mean freeze- 
out time, where T the central temperature at the mean freeze-out time while T s is 
the temperature of the surface of the fireball. Also we introduce a lengthscale, R s , 
defined by the transverse coordinate where the temperature distribution decreases 
to half, as compared to the center, so T(r x = r y = R s ) = To/2. In this notation, a 
flat temperature profile corresponds to the R s — > oo limit. 

Furhermore, rj s is the space-time rapidity of the saddle-point (point of max- 
imal emittivity), which depends on the rapidity, the longitudinal expansion, the 
transverse mass and on the central freeze-out temperature rj s — y = 1 1 mt , where 

y = 0.5 log( ) is the rapidity. Because of this, at midrapidity, where PHENIX 
and STAR V2(j>t) data was taken, Eq. {JJ simplifies to Wit = mt- These definitions 
imply that the effective temperature is transverse mass and rapidity dependent both 
in the reaction plane and out of the reaction plane direction, 

The effective slope parameter, T*, of the asimuthally averaged single particle 
spectra can be defined as 



1 f 1 



T 2 V T T 



(8) 



A relativistic generalization of the kinetic energy, relevant for the elliptic flow 
analysis is defined as 

Finally, the eccentricity of the single particle spectra can be defined as 

(10) 



T X Ty 
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The universal scaling variable w turns out to be the product 

w = ^e (11) 

It is interesting to note, that each of the kinetic energy term, the effective temper- 
ature T* and the momentum space eccentricity e are transverse mass and rapidity 
dependent factors. However, for m t /To ^> 1, e — * e p = ^ 2 ~^ 2 becomes independent 
of transverse mass and rapidity. 

This structure of the universal scaling variable suggests that the detailed cen- 
trality, transverse mass and rapidity dependence of the elliptic flow is only apparent, 
and the Buda-Lund model predicts that when plotted against the scaling variable 
to, a data collapsing behavior sets in and a universal scaling curve emerges, which 
coincides with the ratio of the above two Bessel functions. Additional details about 
the ellipsoidally symmetric model and its result on V2(n) can be found in Ref. [ I15| . 

Eq. (J2J depends, for a given centrality class, on rapidity y and transverse mass 
mt. Before comparing our result to the V2(n) data of PHOBOS, we thus performed 
a saddle point integration in the transverse momentum variable and performed a 
change of variables to the pseudo-rapidity n = 0.5 log( ctz^ 2 ), similarly to Ref. [ 
126) . This way, we have evaluated the single-particle invariant spectra in terms of 
the variables n and 0, and calculated V2(n) from this distribution, a procedure 
corresponding to the PHOBOS measurement described in Ref. [^. 

Transverse momentum dependent elliptic flow data at mid-rapidity can also be 
compared to the Buda-Lund result directly, as it was done in e.g. Ref. 
Comparison to experimental data. When fitting only v^ipi) data of identified 
particles as measured by STAR and PHENIX, we have found that fit parameters u t 
and e are essential. When analysing V2 (77) data, we found that they can be described 
by a two-parameter fit: parameter e controls the peak while parameter An controlls 
the widht of these distributions. Note that the same An parameter controlls the 
width of the pseudo-rapidity distribution as well, but these simultaneous fits will 
be reported elsewhere. The quality of our fits was found to be insensitive to the 
precise value of other parameters of the model. Hence we have fixed their values. We 
also excluded points with large rapidity from lower center of mass energies v-zirj) 
fits (n > 4 for 19.6 GeV, r\ > 4.5 for 62.4 GeV) and points with large transverse 
momentum (p t > 1.3 GeV) from V2(pt) fits. 

Fits to PHOBOS data of Ref. [|Tj, PHENIX data of Ref. [ 2 and STAR data of 
Ref. are shown in the top three panels of Fig.^ Bottom panel of Fig.^demon- 
strates that the investigated PHOBOS, PHENIX and STAR data points follow the 
theoretically predicted scaling law. The fitting package is available at Ref. [ I27|. 
and the fits can be redone on the interactive online Buda-Lund page [I28|. 

In contrast, Ref. [E] considered the description of the centrality and pseudo- 
rapidity dependence of the elliptic flow in 200 GeV Au+Au collisions as measured 
by the PHOBOS collaboration in Ref. [^, however, this attempt failed to describe 
these data in a unified framework and concluded, that "whether all of the observed 
discrepancies between elliptic flow data and . . . simulations can be blamed on 'late 
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hadronic viscosity' . . . depends on presently unknown details of the initial state of 
the matter formed in heavy-ion collisions" . 

This failure can be contrasted to the success of the Buda-Lund hydro model, 
that analytically predicted the scaling properties of elliptic flow in 2003. The par- 
ticle mass, transverse momentum, pseudorapidity and centrality dependence of the 
elliptic flow at RHIC arc being fully explored only now in experimental studies. We 
think that one of the important elements of success was that the Buda-Lund hydro 
model was based on exact, parametric and analytic solutions of hydrodynamics. 
These solutions were found recently for both non-relativistic and relativistic kine- 
matic domain. Non-relativistic solutions include refs. [ [231 EU ESI- Relativistic 
solutions include not only the Hwa-Bjorken and the Hubble solutions but their new 
relativistic generalizations for 1+1 dimensional, 1+3 dimensional axially symmetric, 
and ellipsoidally symmetric solutions, see refs. [E3 I17j . 

Conclusions. In summary, we have shown that the excitation function of the 
transverse momentum and pseudorapidity dependence of the elliptic flow in Au+Au 
collisions is well described by the Buda-Lund model. We have thus provided a quan- 
titative evidence of the validity of the perfect fluid picture of soft particle production 
in Au+Au collisions at RHIC. This perfect fluid extends far away from mid-rapidity. 
The universal scaling of PHOBOS v 2 (rj) and PHENIX and STAR v 2 (pt), expressed 
by Eq. J5J and illustrated by Fig. provides a successful quantitative as well as 
qualitative test for the appearance of a perfect fluid in Au+Au collisions at various 
colliding energies at RHIC. 
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Fig. 1. STAR [0 (top left), PHENIX [0 (top right) and PHOBOS [Q] (mid- 
dle) data on elliptic flow, i>2, plotted versus p t and rj and fitted with Buda-Lund 
model. Bottom: elliptic flow versus variable w is plotted. The data points show the 
predicted [E| universal scaling. 
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PHOBOS Au+Au v 2 (r]) 


19.6 GeV 


62.4 GeV 


130 GeV 


200 GeV 


e 


0.30 ±0.03 


0.37 ±0.01 


0.40 ±0.01 


0.42 ±0.01 


At] 


2.0 ±0.3 


2.28 ±0.05 


2.60 ±0.04 


2.70 ±0.04 


X 2 /NDF 


1.9/11 


21/13 


36/15 


26/15 


conf. level 


100% 


19% 


0.78% 


6.2% 




PHENIX 200GcV Au+Au v 2 (p t ) 


STAR 


TT 


K 


P 


7T, K, p 


e 


0.31 ±0.003 


0.32 ±0.04 


0.3 ±0.1 


0.38 ±0.01 


u t 


1.2 ±0.1 


0.7 ±0.12 


0.6 ±0.3 


1.0 ±0.1 


X 2 /NDF 


21/7 


5/5 


10/2 


45/21 


conf. level 


2% 


80% 


8% 


0.5% 



Table 1. Results of fits to data of PHOBOS [0, PHENIX [ 2 and STAR [0. 
Fixed parameters: To = 200 MeV, a =1, and for V2(r)), Ut = 1.5. Note that data 
from PHOBOS, PHENIX and STAR have different centrality selection. 
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